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Abstract 

Dynamic response of the F = 2 spinor Bose-Einstein condensate (BEC) un- 
der the influence of external magnetic fields is studied. A general formula is 
given for the oscillation period to describe population transfer from the initial 
polar state to other spin states. We show that when the frequency and the 
reduced amplitude of the longitudinal magnetic field are related in a specific 
manner, the population of the initial spin-0 state will be dynamically localized 
during time evolution. The effects of external noise and nonlinear spin ex- 
change interaction on the dynamics of the spinor BEC are studied. We show 
that while the external noise may eventually destroy the Rabi oscillations and 
dynamic spin localization, these coherent phenomena are robust against the 
nonlinear atomic interaction. 
PACS number(s): 03.75.Fi, 05.30.Jp, 32.80.Pj 

Key words: Spinor Bose-Einstein condensate, Rabi oscillation, Dynamic lo- 
calization 



I. INTRODUCTION 



Recent experiments on 23 Na condensates confined in an optical trap [0,0,0 have stim- 
ulated extensive interest in the study of multicomponent spinor Bose-Einstein condensates 
(BECs). Due to the hyperfine spin of atoms, BECs of alkali-metal atoms have internal 
degrees of freedom which are frozen in a magnetic trap. Introduction of an optical trap 
liberates them to allow BEC to be in a superposition of magnetic sublevels. Therefore, 
the spinor BEC can be represented by a vector order parameter. The F — 1 spinor BEC 
was first theoretically studied by Ho [ffl, Ohmi and Machida 0] by generalizing the Gross- 
Pitaevskii equation under the restriction of gauge and spin-rotation symmetry. Within the 
mean-field theory, they predicted a rich set of new phenomena such as spin textures and 
topological excitations. Law et al. § constructed an excellent algebraic representation of 
the F — 1 BEC Hamiltonian to study the exact many-body states, and found that spin- 
exchange interactions cause a set of collective dynamic behavior of BEC. Since the spinor 
BEC appears feasible by using the F = 2 multiplet of bosons such as 23 Na, 87 Rb, or 85 Rb, 
it is necessary to investigate the ground-state structure and magnetic response of F = 2 
spinor BEC. Recently, Ciobanu et al. [|7| generalized the approach for the F — 1 spinor BEC 
developed by Ho to study the ground-state structure of the F = 2 spinor BEC. They found 
that there are three possible phases in zero magnetic field, which are characterized by a pair 
of parameters describing the ferromagnetic order and the formation of singlet pairs. From 
current estimates of scattering lengths they also found that the spinor BEC of 87 Rb and 23 Na 
have a polar ground state, whereas those of 85 Rb and 83 Rb are cyclic and ferromagnetic, 
respectively. Koashi and Ueda || studied the exact eigenspectra and eigenstates of F = 2 
spinor BEC. They found that, compared to F — 1 spinor BEC, the F = 2 spinor BEC 
exhibits an even richer magnetic response due to quantum correlations among three bosons. 

Different from superfluid 3 He, a new feature in a spinor BEC is that its response to an 
external magnetic field is dominated by the electronic rather than nuclear spin. This opens 
up possibilities of manipulating the magnetism of superfluid vapors. Observation of spin 
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domains by Stenger et al. at MIT || offers a remarkable example of such manipulations. Pu 
et al. |J investigated the effects of external magnetic fields on the dynamics of F — 1 spinor 
BEC and found various magnetic field-induced effects, such as stochastization in population 
evolution, metastability in spin composition, and dynamic localization in spin space. In this 
paper we investigate the time evolution of the F = 2 spinor BEC in the presence of external 
magnetic fields with longitudinal and transverse field strengths B z (t) and B x , respectively. 
Here the longitudinal magnetic field lifts the energy degeneracy of the spin states through 
Zeeman effects, whereas the transverse field appears as a coupling between different spin 
components. We assume the longitudinal field to be time-dependent and that the transverse 
field is static. Using perturbation theory, a general formula for the oscillation period to 
describe the passing across the initial spin-polarized state is obtained. Another aspect we 
are interested is the effects of the external noise and the nonlinear spin-exchange interactions 
on the quantum coherent behavior of the spinor BEC. 

In Sec. II the Hamiltonian for the spin-2 BEC system is presented. The Rabi oscillation 
and dynamic spin localization is shown in Sec. III. In Sec. IV we discuss the effects of 
external noise on the dynamics of the spinor BEC. A full numerical discussion based on an 
effective one- dimensional nonlinear Schrodinger equation is given in Sec. V. A summary is 
given in Sec. VI. 

II. THE MODEL: SINGLE-MODE APPROXIMATION 

We consider the F = 2 spinor BEC subject to a spatial-uniform magnetic field B(t). 
The Hamiltonian invariant under spin space rotation and gauge transformation is written 
in terms of the five- component field operators: \l/ + 2, ^-2, corresponding to the sublevels 
rrip = +2, —2 of the hyperfine state F = 2. Namely, it is given |7||| by H = H + H B , 

Ho = J <M^W+ • V*. + U^ a + +| A) 2 

i 

+c 2 *+#+(2a;2a' | 00) x (00 | 2(3; 2/3') W], (1) 
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(2) 



where c , 5i, and c 2 are related to scattering lengths a , a 2 , and a 4 of the two colliding bosons, 
with total angular momenta 0, 2, and 4, by Co = 47r/i 2 (3a4+4a2)/7M, ci = knf^ia^ — a^jlM^ 
and C2 = 4:7iH 2 (3a4 — 10a2 + 7a )/7M. In addition, we have introduced in Eqs. (1) and (2) 
5x5 spin matrices Fj (i = x, y, z) 
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We assume the external magnetic field is weak so that the coordinate wave function $(r) 
is independent of the spin state and solely determined by the first three terms of Eq. (1), 
namely, [-h 2 V 2 /2M + U + c (iV- 1) |$| 2 ]$ = with iV being the total particle number. 
Substituting = a a $ into Eqs. (1) and (2) and keeping only spin-dependent terms, we 
obtain 



and 

H B = -fi, B g f B(t)-F, (7) 

where Fi = a^(Fi) a pap, S + = = (Eg ) 2 /2 — afati + a^a^, and q = (q) / |$| 4 . We 
further assume that the magnetic fields consist of longitudinal and transverse components. 
Without loss of generality, the transverse direction of the field is chosen to be along the x 
axis, i.e., B{t) = Bi(t)z + B x x. In such a case, the second-quantized Hamiltonian of H B is 

Hb = —HB9fBi(t)(a,2a2 + afai - cT^S-i - a± 2 a_ 2 ) 

—fjLBgfB x (a,2ai + ^3/2afa + y3/2a,Q a_i + al 1 a_ 2 + h.c). (8) 

Similar to Ciobanu et al. , mean-field approximation is used such that the field operators 
a a are replaced by c numbers a a = \/P^e lda , where P a = N a /N is the population in spin 
a and 9 a the phase of wave function a a . Furthermore, since this paper deals with the 
quantum coherent behavior of the system under the influence of the external magnetic 
fields, we assume that the initial spin state of the BEC is the eigenstate in the absence of 
external fields, consequently, contribution from Hamiltonian (6) is a constant energy shift, 
and can be neglected in the dynamics. The semiclassical equations of motion in Heisenberg 
representation can be derived from the Hamiltonian Hb (We introduce the state vector 
a = (a 2 , ...,a_ 2 ) T ) 

i a= H e ff(t)a, (9) 

where 

H eff {t) = -k{t)F z - b x F x . (10) 

Here we have defined bi(t) = ^B9fBi{t) and b x = ^B9fB x . The time evolution of the system 
has a rich set of intriguing features with different initial conditions. In this paper, we consider 



the case that the system begins with unperturbed spin-0 state a(0) = (0, 0, 1, 0, 0) T , i.e., the 
P0 polar phase defined by Ciobanu et al. J7|. This phase is possible to be realized in 87 Rb 
and 23 Na. 

III. RABI OSCILLATIONS AND DYNAMIC SPIN LOCALIZATION 

In this section, we analyze the time evolution of F = 2 BEC described by Eq. (9) with 
initial polar state a(0) = (0, 0, 1, 0, 0) T . Assuming that a(t) = e*Jo dThli - T)Fz ip(t) , we obtain 

i<p(t) = Hj(tMt), (11) 

where Hi(t) = — 6 x exp(— i J*J drbi(r)F z )F x exp(i J * drbi(r)F z ). Making use of the identity 

exp(i\F z )F x exp(—i\F z ) = F x cos A — F y sin A, (12) 

the formal solution of Eq. (11) is obtained 

<p(t) = fe'fo ^[X(T)F x +Y(r)F y](p ( ^ (13) 

where X(t) = b x cos[/q drbi(r)], Y(t) = b x sin[/o drbi(r)], and T denotes time ordering. Con- 
sidering the situation in which the transverse magnetic field b x is weak, we may approximate 
Eq. (13) by the following solution 

<p(t) = jfi Mx(T)F x +nT)F y](p ( y (M) 

This solution is valid to the first order in b x and preserves unitarity. After a straightforward 
calculation, we obtain the time evolution of the population of spin-a state as 

P a (t) = \a a (t)\ 2 = 

where A = \J (J * drX{r)) 2 + (/ * <irF(r)) 2 , and the matrix element is defined as 

d 2 pa (9) = (20 | e-*^ | 2a), (P, a = 2, -2). (16) 
In particular, the population of the initial spin-0 state at time t is 
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P {t) = \a (t)\ 2 = cos 2 



JO 



(17) 



Thus to the first-order approximation in b x we obtain the analytical expression of the evo- 
lution probability for the spinor BEC remaining in the initial polar state. Expression (17) 

, "'' T ' can be 

expanded as a discrete Fourier series, and the time integral is either bounded or increases 
linearly in time (on top of an oscillatory piece). In the former case, Pq remains close to 
unity at all times because of the smallness of b x . This high population of the initial spin 
state is just the spin localization. In the latter case, P oscillates between and 1, imply- 
ing population transfer between P and other P a (a = ±1, ±2). Hence we have the spin 
derealization, i.e., Rabi oscillation. Thus, we obtain the oscillation period T as 



T = 



7T 



lim 



f t dr 1 b x e i fo ldT2bl{T2) 



(18) 



This result is still valid for any time-dependent external fields. We give two examples of 
special time-dependent external magnetic fields. 

(a) Sinusoidal field bi(t) = bcos(ujt). According to Eq. (18), the oscillation period is 



b x J Q {b/uY 

where J is the zeroth order Bessel function of the first kind. From this result, it can be seen 
that when J (b/u) ^ 0, there exists a finite oscillation period, indicating that the population 
can transit from an initial polar state to other spin states within the driving process. When 
6 = 0, the longitudinal field-free oscillation has period T = n/b x . The fact \J Q (b/uS)\ < 1 
implies that T > T . This shows that the invasion of the longitudinal magnetic field will 
suppress the coherent population transfer among the five spin components. The extreme 
case occurs when the reduced amplitude of the longitudinal magnetic field &/w is a root 
of Jo- In this case, we have T — > oo, which implies that the population transfer is totally 
suppressed. Hence the system will stay in the initial polar state during the whole driving 
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process. This is just the phenomenon of spin localization. The above analysis is illustrated 
in Fig. 1, where the scaled oscillation period T/T is plotted against b/u. It can be seen 
that the oscillation period increases significantly when b/u approaches the roots of Bessel 
function J . 

(b) A combination of static and sinusoidal magnetic fields bi(t) = b + bcos(ut). In this 
case, the oscillation period is 

T = . J „, ; ^ , (20) 



b x lim^co SSS 4/,(^) 



where J is the Anger function |10[ defined by 



— 1 f n 

J a (b) = — / dx cos(ax — b sin x) . (21) 

7T JO 

There exists an infinite oscillation period when bo/u ^ k (k is an integer), which means 
that the population time evolution will undergo localization in the initially populated spin-0 
state during the whole driving process. When b /u = k, the oscillation period becomes 

T = b x J k (b/u)' (22) 

where are k-th order Bessel functions. From this result, we find that we can have a finite 
oscillation period only when b/u is not equal to a root of J^, and in this case the evolution 
of the initial polar phase is delocalized into other spin states. If b/u becomes a root of J^, 
the oscillation period approaches infinity, and we have the spin localization. 



IV. EFFECTS OF EXTERNAL NOISE 

In practice, since the external magnetic field may have a fluctuation component, the 
effects of external noise has to be considered. Another motivation for the introduction of 
an external noise is to provide the fluctuations required to destroy the coherence of the 
population transfer process and, in suitable cases, leads to a rate process for the decay of 
the population. In the following, we study the dynamics of the system when the external 
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magnetic field has a fluctuating component. In the simplest case, the statistical properties 
of the imposed magnetic field can be prescribed to be independent of the characteristics of 



the system. As a versatile choice for the noise, an Ornstein-Uhlenbeck (OU) process [II] is 
used. This permits the investigation of the role of the strength and size of correlation time 
of the noise on the time evolution of the spinor BEC. Without loss of generality, the external 
longitudinal magnetic field bi(t) is assumed to contain two components, a stochastic part 
fit) and a systematic part 60 + &i cos(uit), i.e., 

hit) = b + b lC os(ujt) + f(t). (23) 

The noise f(t) is assumed to be characterized by an OU process whereby it has zero average 
value and correlation function 

(f(t)f(s))=A 2 exp(-\t-s\/r c ). (24) 

Here the quantities A and r c are the strength and decay constant of the noise, respectively. 
When the noise is external to the system, and therefore not necessarily thermal in character, 
A and r c can be varied in a controlled manner, and are not restricted by the physical 
properties of the system. In the absence of an analytic solution, a numerical investigation 
of population P a (t) is adapted. The population distribution in the five components of 
spinor BEC with the external noise are solved numerically by generating trajectories for the 
different realization of the noise. The procedure |L2| to integrate the stochastic equations is 



as follows. A stochastic term is added at each step with its statistical properties described by 
an OU process. The OU process is generated by solving a Langevin equation with a delta- 
correlated noise term. This ensures that the correlation function of f(t) has the desired 
statistical property given by Eq. (24). 

We present in Fig. 2 the time evolution of the population Po(t) of spin-0 state in the static 
magnetic field \b\ = in Eq. (23)] for the value of bo = 10b x . In the absence of noise, as shown 
in Fig. 2(a), the population of spin-0 state is always close to unity during time evolution. 
This effect is purely the static spin localization. In this case, the longitudinal magnetic field 
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lifts the degeneracy of the five spin components and consequent energy mismatch hinders 
the BEC transfer from the initially occupied polar state to its neighbours in the spin space. 
In the case of the intermediate (At c ~ 1) noise modulation, as shown in Fig. 2(b) (solid 
curve), it can be seen that the population of spin-0 state has a damped oscillation with a 
slow monotonic decay superimposed. The combination of noise and static magnetic field, 
suppressing population transfer, produces the decay. It also shows in Fig. 2(b) that with 
increasing strength of noise, the decay becomes more rapid. Therefore, the decaying time 
can be controlled by varying the magnetic field or the strength of noise. 

Figure 3 shows the time dependencies of Po(t) in the time- dependent longitudinal mag- 
netic field with the systematic part bi = b cos (ujt). The field parameters are chosen as b = 2b x 
and 00 = 106^, corresponding to the Rabi oscillation with scaled period T = T /J (0.2), as 
shown in Fig. 3(a). In the case of a weak noise A/b x = 0.2 (dashed curve), the population 
Po is still oscillatory in our scope of time, but its oscillation amplitude decreases with time. 
When the strength of the noise increases to A/b x = 1.0 (dotted curve), the Rabi oscillation 
breaks down completely and the system decays fast towards equilibrium. This suggests that 
the coherent Rabi oscillation for a spinor BEC is sensitive to the external noise and is de- 
stroyed even in a weak coupling regime. We also present in Fig. 4 the time evolution of P 
for the value of b/u — 2.402, corresponding to the dynamic localization, as shown in Fig. 
4(a). In Fig. 4(b), we can see that in the case of intermediate (Ar c ~ 1) noise modulation 
(solid curve), the population P still remains close to unity for an extremely long time. This 
insensitivity to the presence of weak noise reflects the strength of the systematic field. When 
the strength of noise increases to a strong coupling regime, a more rapid and less oscillatory 
decay is observed in the population evolution, as shown in Fig. 4(b) (dotted curve). 

Figure 5 shows the time evolution of Po in the external longitudinal magnetic fields 
consisting of static and time-periodic components [Eq. (23)]. The field parameters are chosen 
as bo/b x = 10, ujjb x = 10, and b\jb x = 38.3, corresponding to the spin dynamic localization, 
as shown in Fig. 5(a). In the case of intermediate (Ar c ~ 1) noise modulation, similar to 
that shown in Fig. 4(b), the population P of spin-0 state remains close to unity in our 
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scope of time. Even in the strong coupling regime (dashed curve), the fast oscillations still 
exist, suggesting the fundamental interplay between the systematic and noise fields. With 
further increase of the strength of noise, the system decays rapidly towards equilibrium, and 
therefore, the spinor BEC initially localized in spin-0 state is delocalized and diffuses among 
all spin states, as shown in Fig. 5(b) (dotted curve). 

V. DYNAMIC SPIN LOCALIZATION BEYOND SINGLE-MODE 

APPROXIMATION 

In the above discussions, the single-mode assumption (SMA) has been used: the wave 
function for each spin component retains the same spatial profile during the time develop- 
ment. Therefore, the possibility of high-energy mode excitation is neglected and the system 
dynamics is fully described by the internal population transfer among the spin components. 
This SMA treatment is valid for a small particle number N. When iV is large, the nonlinear 
spin-mixing processes induce the large energy transfer from the initial ground state to the 
excited modes, thus destroying the validity of the SMA. However, when the initial state is 
taken to be the eigenstate of the mean-field G-P equation, the SMA is expected to describe 
well the dynamics in the absence of dissipation and external driving |13| . In this section, we 
go beyond the SMA by means of a numerical simulation of the time-spatial evolution of the 
system. 

Due to anisotropic nature of the optical trapping potential, the cigar shaped BEC is 
assumed to be quasi-one-dimensional and hence the wave function associated with the spin- 
et state may be written as 

i/> a (x, y, z, t) = <t> ± (x, y)(j> a (z, t)e~^\ (25) 

where z is the direction of weak confinement, <p±(x,y) is the ground-state wave function of 
the two-dimensional harmonic oscillator, and uj± the tight confinement frequency. Inserting 
Eq. (25) into Eqs. (1) and (2), we obtain the equations of motion for the longitudinal wave 
functions <f) a (z,t) in a dimensionless form 
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ih^ = £^2 + f-ipi + fzi>2 + ^* 2 s_ - Mi - 26,(0^2 (26a) 



ifi^ = Ztyi + 7+^2 + sl\f-^ + fz^i - r-is- - 6x(^i + t/^o) - 6i(*)^i (26b) 



ifi^ = + + f-*f>-i) + r o s- - b x (^ + (26c) 



ifi^=± = Zty-i + + f-ip-2 - M-i - #fs_ - My + ^-2) + 6i(*)^-i (26d) 

= £^ 2 + - fA-2 + r 2 s- - M-i + 26,(0^-2 (26e) 

where £ = -rf 2 /rfz 2 + 2; 2 /4 + CoiV?7(|02| 2 + |0i| 2 + |0o| 2 + |0-i| 2 + 10-2| 2 ) with transverse struc- 



ture factory = / \(p ± ( y x,y)\ A dxdy/ j {(p^x^^dxdy, /+ = c 1 (^0i+y / 3/20*0 o +y / 3/ 2 0( l S0-i+ 
0*^-2), /- = /;, / 2 = (ci/2) E-=^ 2 |0,| 2 , and s_ = (2c 2 /5)(0 2 /2 - fafa + 2 0_ 2 ). The 
above equations have been written in dimensionless form and the units for length, energy, and 



time are ^h/(2muj z ), fku z , and l/u z , respectively, where oj z is the axial trapping frequency. 
In the numerical simulation of the dynamic equations (26), the initial wave functions are 
taken to be ground state of the one-dimension G-P equations C(f>i(z,0) = fj,(f>i(z,0), where 
(f>i(z, 0) satisfies the normalization / \<f>i(z, 0)\ 2 dz = Ni/N with iVj the initial particle number 
in spin-i component. 

We first discuss the dynamics of the system in the absence of coupling magnetic fields. 
Figure 6 shows the time evolution of the population of spin-0 component for two kinds of 
initial population. The solid line in Fig. 6 corresponds to the initial state that five spin 
components are initially equally populated, whereas the dotted line corresponds to the case 
that the BEC initially populates spin-0 components. It shows (solid line) that during time 
evolution, the population of spin-0 component develops into a complex oscillatory structure, 
implying the excitation of high-energy modes. The excitation originates from the fact that 
the initial state is not a mean- field eigenstate of the SMA Hamiltonian (6), thus leading 
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to a redistribution of the total energy E between the symmetric and unsymmetrical parts 
(denoted by E s and E a , respectively) of the total Hamiltonian. When E a is sufficiently 
large, the higher modes are excited by the nonlinear spin-mixing processes and the SMA 
description breaks down for large value of particle number N. When the initial state is 
chosen to be the SMA eigenstate, the population of spin-0 state does not vary in time 
(dotted line). Thus, SMA treatment is appropriate for the spin polarized state with the 
BEC localized in spin-0 component. 

In the presence of external magnetic fields, the five spin components are coupled. The 
atomic population is featured by Rabi oscillations among the spin states. To see the validity 
of the SMA in studying the dynamic spin localization, We present in Fig. 7 the time evolution 
of the population of spin-0 state for several values of amplitude of oscillatory magnetic 
field. It shows that the weak transverse magnetic field couples the five spin components, 
inducing Rabi oscillation between the initial populated spin-0 state and the other spin states. 
The presence of an additional oscillatory magnetic field slows down the oscillation periods. 
Finally, when the field parameter is modulated to satisfy b/ftu = 2.4, as shown in Fig. 7(b), 
the Rabi oscillation completely stops and population of spin-0 state always remains close 
to unity during its time development. The good agreement between the SMA results and 
exact numerical solution suggests that the phenomenon of dynamic spin localization is well- 
defined beyond the simple SMA. This robustness of dynamic localization against nonlinear 
atomic interaction is contrast to the case of two-component BEC system, where there is no 
energy-stable spin-polarized states. 

VI. CONCLUSION 

In summary, by the application of the external magnetic fields which may contain a 
random component, we have elaborated on the control of population transfer among the 
internal spin states for the F = 2 spinor BEC. After taking the mean-field approximation a 
general formula is obtained for the oscillation period to describe the quantum transition from 
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the initial spin polar state to other spin states. In particular, when the frequency and reduced 
amplitude of the longitudinal magnetic field are related in a specific manner, the population 
of the initial spin-0 state is dynamically localized during the time evolution. In addition, the 
effects of external noise on the coherent population transfer dynamics of the spinor BEC are 
studied. We found that (i) The Rabi oscillations are sensitive to the external noise, and even 
weak coupling can destroy completely the Rabi oscillations in an earlier length of time, (ii) 
When the strength of noise is weak, because of strong suppression caused by the systematic 
field, the static and dynamic localization may remain for experimentally relevant times. 

(iii) When the strength of noise is strong, derealization due to phase randomness occurs 
during time evolution and the initially localized BEC diffuses among five spin components. 

(iv) The Rabi oscillation and dynamic spin localization is robust against the nonlinear spin 
exchange interaction. Thus in order to observe the static or dynamic localization in spin 
space one should produce stable magnetic fields, which has as little as possible the stochastic 
component. Although it is difficult to keep the external noise from intruding on dynamics, 



recent observation of the well-defined coherence between the two BECs suggests that 

the external noise can be greatly suppressed in experimental setup. Therefore, we expect 
such macroscopic coherent quantum controls as Rabi oscillation and dynamic localization 
may be realized in spinor BEC experiments. 
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Figure captions 

Fig. 1. Scaled oscilltion period T/T of the sinusoidal magnetic field case, versus the 
reduced field strength b/uo. 

Fig. 2. Time evolution of Po(t) with and without external noise in the static magnetic 
field, b/b x = 10. (a) the case without the noise, (b) the case with the noise A/b x = 1.0 (solid 
curve), A/b x = 5.0 (dashed curve), and A/b x = 10.0 (dotted curve). b x r c is fixed to be 1.0 
for the case with the noise. 

Fig. 3. Time evolution of Po(t) with and without external noise in sinusoidal magnetic 
field, b/b x = 2, uj/b x = 10. (a) the case without the noise, (b) the case with the noise 
A/b x = 0.2 (solid curve), A/b x = 0.6 (dashed curve), and A/b x = 1.0 (dotted curve). b x r c 
is fixed to be 1.0 for the case with the noise. 

Fig. 4. Time evolution of Po(t) with and without external noise in sinusoidal magnetic 
field, b/b x = 24.02, uo/b x = 10. (a) the case without the noise, (b) the case with the noise 
A/b x = 1.0 (solid curve), A/b x = 5.0 (dashed curve), and A/b x = 10.0 (dotted curve). b x r c 
is fixed to be 1.0 for the case with the noise. 

Fig. 5. Time evolution of Po(t) with and without external noise in a combination of 
static and sinusoidal magnetic field, bo/b x = 10, b\jb x = 38.3, uj/b x = 10. (a) the case 
without the noise, (b) the case with the noise A/b x = 1.0 (solid curve), A/b x = 5.0 (dashed 
curve), and A/b x = 10.0 (dotted curve). b x r c is fixed to be 1.0 for the case with the noise. 

Fig. 6. time evolution of Po(t) for two kinds of initial population: Ni/N = 0.2, i = 
—2,..., 2 (solid line), N /N = 1 (dashed line). Other parameters are oo z = 2ir x 60Hz, 
a = 34.9a,e, a 2 = 45.8a_B, a 4 = 64.5a_B for 23 Na (as is the Bohr radius), rj — 1, N — 10000. 

Fig. 7. Time evolution of Po(t) in the presence of a uniform field b x = 1. (a) the 
case without sinusoidal magnetic field, (b) the case with sinusoidal magnetic field bi(t) = 
bcos(ut), uj = 10, b = 24. Other parameters are the same as used in Fig. 6. 
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